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FAST NEUTRON RESONANCE RADIOGRAPHY FOR ELEMENTAL MAPPING 

GOVERNMENT SUPPORT 

The invention was supported, in whole or in part, by the Office of National Drug 
Control Policy (ONDCP), under Contract DAAD07-98-C-01 17, and by the Federal Aviation 
5 Administration, under Grant 93-G-053. The US Government has certain rights in the 
invention. 

BACKGROUND 

Today, both narcotics trafficking and terrorist attacks are worldwide problems. Drug 
dependence is a chronic, relapsing disorder that exacts an enormous cost on individuals, 
10 families, businesses, communities, and nations. Terrorist highjacking and bombings, 

especially those involving and targeting civilian airplanes, greatly jeopardize the lives of the 
public. 

Drugs, explosives and weapons are transported illegally by various methods. Large 
objects such as trucks or cargo containers are checked at the United States' southwest border 
1 5 with Mexico for large amounts of drugs or other goods of interest, while passenger's luggage 
is screened for a small amount of explosives or drugs (1 00—300 grams) or for weapons at 
airports. There are a number of techniques currently used or proposed for defecting drugs, 
explosives and/or other contraband. 

There have been many governmental efforts in the United States to develop detection 
20 technology for drugs, explosives and weapons led by the Office of National Drug Control 
Policy (ONDCP) and the Federal Aviation Administration (FAA). Commonly used methods 
for detecting drugs, explosives and/or weapons in parcels or cargo are described, below. 

The oldest and perhaps the most reliable way of detecting hidden drugs and explosives 
is to check the suspected object manually. After the baggage or cargo container is opened and 
25 carefully searched, we can determine with a high level of confidence whether a weapon, an 
explosive device or an illegal drug is present. However, manual checking is slow, expensive 
and perhaps dangerous. In addition, manual inspection can violate societal norms and legal 
standards relating to privacy. Usually, only some randomly selected or highly suspected items 
are checked manually. 

30 Use of a detector dog and its handler is another old, but still broadly used, technique 

for drugs and explosives detection. Canine olfactory systems are capable of reacting to the 

1 
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substance, it learns to discriminate the vapor of that substance from other odors in the 
environment by reacting to the compound(s) and earning positive reinforcement from the 
handler when it does so. A dog usually works for several hours a day but may not be willing 
to work at the required time. In addition, canine detection fails if no vapor is coming out of 
5 the package. 

Drugs and explosives can also be detected by trace contaminants such as vapor or 
particulate residue. Air samples are collected, condensed, and then analyzed, very often with 
ion mobility spectrometer (IMS) detectors. Ways of collecting samples range from sniffing 
the vapors to wiping the surface, but they all depend on contamination. This detection method 
10 is used to screen people and/or boarding passes. 

X-ray detection systems are the most commonly used for parcel and cargo inspection. 
X-ray transmission radiography or tomography is based on the attenuation of x-rays to image 
the electron density in the object under inspection. An experienced operator can identify 
different items by their density and shape distribution. X-ray systems are very effective in 
1 5 telling metal from organic material, but lack the ability to distinguish drugs or plastic 

explosives from ordinary goods. When low energy x-rays are used as the source, a second 
detector may be used to record the scattered x-rays. This generally helps to map low-Z 
material near the surface of the inspected object and is used to detect sheet explosives hidden 
in the parcel wall or other contraband beneath the wall of a cargo container. To reduce the 
20 cost and increase penetration, a radioactive gamma-ray source ( l37 Cs or ^Co) is sometimes 
used to replace the x-ray machine. 

The use of neutron resonance radiography (NRR) has been described for applications 
other than detection of items such as explosives, weapons and drugs. NRR is based on the fact 
that, for a given element, the absorption of a beam of neutrons varies with neutron energy. 
25 Neutron absorption for a given element will typically show a series of relatively narrow peaks 
("resonances' 1 ) as a function of energy. If an image is taken at an energy matching the 
resonance energy and another image is taken at an energy off the resonance energy for a 
particular element, the difference of the two images will provide an enhanced image of that 
element. Thus, with appropriate choices of energies, an image may be taken of a particular 
30 element (e.g„ carbon) within the object 

Neutrons of appropriate energies may be produced using a particle accelerator to 
accelerate beams of charged particles, such as protons or deuterons, onto a solid or gas target. 
In the usual method of performing resonance radiography, the energy of the accelerator is 
varied so as to produce beams of neutrons with a narrow range of energies. Each element is 
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imaged at two energies. Further, the target used is made relatively thin in order to produce 
neutrons with a narrow range of energies. Such thin targets are inefficient producers of 
neutrons and also limit the accelerator output due to heating in the target. 

SUMMARY 

5 Methods and apparatus for mapping the elemental composition of objects, particularly 

parcel-size objects, by means of neutron resonance radiography (NRR) are further described, 
below. In methods of this disclosure, neutrons are produced in a relatively broad band of 
energies from an efficient thick target, and the angular dependence of the neutron energy can 
be used to scan multiple energies simultaneously and to measure several elements 

1 0 simultaneously. Further, the energy of the particle accelerator can be fixed; and thus, 
complications involved with changing accelerator energies are avoided. Compact fixed- 
energy accelerators may thus be deployed. 

Unlike other resonance techniques, such as pulsed fast neutron transmission 
spectroscopy (PFNTS), the methods and apparatus disclosed herein do not require an 

1 5 accelerator that produces short pulses of neutrons; moreover, these same methods and 

apparatus can use a simple large imager as a detector. Complexity in both the accelerator and 
the detector is consequently avoided. In contrast, the PFNTS method determines neutron 
energy by measuring the velocity of the neutrons, and thus requires both a very short (on the 
order of nanoseconds) pulse of neutrons and a relatively long distance over which to measure 

20 their time of flight so as to determine energy, thereby increasing the space required for such a 
system. 

In NRR methods further described below, two-dimensional elemental mapping of, for 
example, hydrogen, carbon, nitrogen, oxygen and the sum of other elements can be calculated 
using fast neutron radiographic images taken at different neutron energies chosen to cover the 

25 resonance features of one or more elements. A radiographic image provides the two- 
dimensional mapping of the sum of elemental contents weighted by their attenuation 
coefficients, and images taken at different neutron energies (at which the elements will have 
different attenuation coefficients) form a set of linear equations for each pixel in the image. 
These equations can then be solved to map individual elemental content. Contraband, such as 

30 explosives and drugs, can be identified by their characteristic elemental composition. 

Accordingly, pieces of passenger luggage can be imaged, in an airport, for example, with the 
methods taught herein to detect contraband such as explosive materials or drugs (e.g., cocaine) 
by their characteristic elemental compositions. Neutrons having different energies (e.g., 2-6 
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MeV) can be obtained at different angles from a deuteron-deuteron (DD) neutron source. A 
fixed-energy radio frequency quadrupole (RFQ) accelerator with a thick target can be used as 
the neutron source in NRR. 

BRIEF DESCRIPTION OF THE DRAWINGS 

5 The foregoing and other features and advantages of the invention will be apparent from 

the following, more-particular description. In the accompanying drawings, like reference 
characters refer to the same or similar parts throughout the different views. The drawings are 
not necessarily to scale, emphasis instead being placed upon illustrating particular principles, 
discussed below. 

10 FIG. 1 charts the neutron cross-section (measured in bams) for hydrogen (H), carbon 

(C), nitrogen (N) and oxygen (0) over an energy range of 0.5 to 2.5 MeV. 

FIG. 2A illustrates experimental transmission of four compositions for neutrons at 2.3 
MeV, wherein different shadings, which are more easily distinguishable in color, represent 
different degrees of neutron transmission and attenuation. 
1 5 FIG. 2B illustrates a ratio, of off-resonance versus on-resonance neutron permeability 

for nitrogen in the four compositions of FIG. 2 A. 

FIG. 2C illustrates a ratio of off-resonance versus on-resonance neutron permeability 
for oxygen in the four compositions of FIG. 2 A. 

FIG. 2D illustrates a ratio of off-resonance versus on-resonance neutron permeability 
20 for carbon in the four compositions of FIG. 2A. 

FIG. 3 illustrates the attenuation, or neutron cross-sections, for hydrogen, carbon, 
nitrogen and oxygen over a range of neutron energies from 2 to 6 MeV. The illustrated boxes 
above the charted attenuation curves represent "energy bins" showing the range of neutron 
energies in this example at various angles shown in FIG. 6. 
25 FIG. 4 illustrates the vector analysis performed using a Householder transform. 

FIG. 5 is an illustration of a fast neutron radiography apparatus. 
FIG. 6 illustrates an apparatus including a deuteron accelerator and a neutron-emitting 
target; the apparatus also includes an object and a detector aligned relative to the target and 
displaceable along an arc about the target to receive neutrons emitted at different angles. 
30 FIG. 7 shows a neutron radiographic simulated image of a bag at 0°, as shown in FIG. 

6, and marks the locations of hidden explosives, weapons and drugs. 

FIG. 8A provides a mapping of the projected hydrogen content in the simulated bag of 

FIG. 7. 
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FIG. 8B provides a mapping of the projected carbon content in the simulated bag of 

FIG. 7. 

FIG. 8C provides a mapping of the projected nitrogen content in the simulated bag of 

FIG. 7. 

FIG. 8D provides a mapping of the projected oxygen content in the simulated bag of 

FIG. 7. 

FIG. 8E provides a mapping of the summed projected content of all other elements in 
the simulated bag of FIG. 7. 

FIG. 8F provides a mapping of the projected summed content of all elements in the 
simulated bag of FIG. 7. 

DETAILED DESCRIPTION 
1. Introduction 

Fast Neutron Resonance Radiography (NRR) is a nuclear technique for the 
determination of two-dimensional elemental content distribution. Offering the major 
advantages of nuclear techniques, NRR can provide elemental information for detection of, 
e.g., explosives, weapons or drugs, and has the ability to penetrate thick objects. 

An accelerator is used to produce source neutrons, which can have energies in the 
range of 2-6 MeV. The source neutrons are directed through an object under examination and 
are detected by a hydrogen-rich scintillator, which typically is formed of plastic. light 
emitted from the scintillator is recorded using a solid-state charged coupled device (CCD) 
camera to form a radiographic image. For each image, the energy spectrum of the source 
neutrons is chosen to cover resonance features of one or more elements of interest. 

The simplest case of NRR is to map one element at a time. FIG. 1 charts the neutron 
cross section (measured in bams) for hydrogen (H) 1 1, carbon (C) 13, nitrogen (N) 15 and 
oxygen (0) 17 over an energy range of 0.5 to 2.5 MeV. A resonance peak/valley (marking a 
high or low point on the vertical axis) for one element is targeted where the cross sections of 
other elements are relatively flat over the same energy range. A radiographic image is taken at 
the "on-resonance" neutron energy the energy at the peak or valley), and another 
radiographic image is taken at a second energy, typically an "off-resonance" neutron energy 
(i.e, an energy level that is neither a peak nor a valley). The pixel-by-pixel difference of the 
two images provides a two-dimensional map of the corresponding element. 
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Several representative materials were tested, and the images produced are shown in 
FIGS. 2A-2D, which illustrate "typical" contraband. The materials were packaged in 3-cm- 
diameter plastic containers, each containing 10 to 30g of material. The circle in the upper left 
of each image represents a drug simulant 12 with the following elemental composition: niH = 
1.5g, mc = 16.3g, m N = 0.9g and m 0 = 5.2g. The circle in the upper right of each image 
represents graphite 14 comprising 10.8g carbon. The circle at lower right represents melamine 
16 comprising the following elements: m H = 1 .8g, mc = 10.6g and m N = 24.9g. Finally, the 
lower-left circle represents an explosive simulant 18 having the following elemental 
composition: m H = l.Og, mc - 7.1g, m N = 7.1g and mo = 12.3g. 

FIG. 2A illustrates experimental transmission of the four compositions for neutrons at 
2.3 MeV. The original color images reflect regions where neutron transmission is greatest 
with the red end of the visible color spectrum, while regions where neutron transmission is 
lowest (and hence, attenuation is highest) are reflected with the blue end of the color spectrum. 
Of the four compositions, melamine 16, the image of which was a blue-tinged green, had the 
highest neutron attenuation in this example; meanwhile, graphite 14, the image of which was 
an orange-red, had the highest neutron transmission. FIGS. 2B, 2C and 2D show respective 
ratios for off-resonance versus on-resonance neutron transmission for nitrogen, oxygen and 
carbon; in the original color images, yellow and, particularly, red areas highlight regions 
where the respective element is detected. 

It becomes increasingly difficult to find a single peak when additional elements are 
considered. Resonance peaks found as described above are generally very narrow; 
accordingly, neutron beams with energies isolated at these peaks are difficult to obtain, even 
when using a variable-energy charged-particle accelerator and a thin target to generate a 
monoenergetic neutron source in accordance with other methods. Moreover, the use of thin 
lithium (Li) targets limits beam currents due to target heating. Further still, neutron energy 
does not drop uniformly across different angles. When objects of finite size are exposed to the 
neutrons, different parts of the object are exposed to neutrons at different angles. 
Consequently, one part of an image may be "on-resonance" while another part of the image is 
not. 

Another approach is to consider each radiographic image as a two-dimensional map of 
projected attenuation. It can also be thought of as a two-dimensional map of the sum of 
content of all existing elements, weighted by their attenuation coefficients. For each pixel in 
the image, there exists a linear equation stating that the total attenuation equals the weighted 
sum of projected elemental contents. When another radiographic image is taken at a different 

6 
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energy spectrum, the resulting linear equation has different attenuation coefficients 
weighting factors) and total attenuation, but has the same projected contents because the object 
is the same. In principle, when there are more equations than the number of existing elements, 
the set of linear equations can be solved for a definite least-squares solution of projected 
5 elemental contents. 

The attenuation coefficients are related to the total neutron cross sections over the 
energy spectrum of the source neutrons. Unlike x-ray attenuation, fast neutron attenuation 
coefficients vary only moderately from energy to energy and from element to element The 
difference is typically within a factor of 2 to 3, which is considered to be neither too small nor 

10 too large to facilitate resolution of the linear equation set. Elements other than hydrogen, 
carbon, nitrogen, and oxygen are modeled as one component to reduce the number of 
equations required. H, C, N, O, and "other" are used to represent these five components. This 
approach is generally appropriate to luggage scanning, though other elements or combinations 
of elements can be mapped, as desired, in accordance with these methods. 

15 By using the broad resonance and linear equation methods described herein, image 

pixels can be calibrated independently, thereby eliminating the need to have the same incident- 
neutron energy across the object. 

2. Fast Neutron Resonance Radiography 

When a radiographic image is taken with neutrons having a broad energy range, (e.g. t 
20 4.487-5.525 MeV, as shown by the 25° stripe in FIG. 3) resonance features of more than one 
element are covered. It turns out that two such images usually are not enough to give us the 
two-dimensional mapping of any element. However, when more than two images are taken, 
the two-dimensional mapping of all elements involved can be solved. 

Using a broad-spectrum neutron source exploits the characteristically broad resonance 
25 features of elements, thereby improving resonance contrast in comparison to choosing narrow 
peaks. In addition, differences among the images (e.g., seven or more) are summed, which 
facilitates mapping of elements (e.g. 9 nitrogen) that are present at low concentrations. 

3. Linear Attenuation Model 

An energy bin that contains a resonance feature of one or more elements is selected, 

30 and a radiograph is taken using neutrons having energies within that energy bin. Each element 

attenuates neutrons as if this single element was alone in the path, and the total attenuation, b, 

is the sum of attenuation by all elements present. The so-formed equation, stating that the 

7 
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total attenuation is the sum of attenuation by all elements, is not sufficient to tell us the 
contents of any element. We thus take radiographic images at additional energy bins and 
obtain additional similar equations. For example, we can take m radiographic image at m 
energy bins for the bulk material composed of n elements, and we obtain the following 
5 equations for each pixel over the image: 

an xj + ai 2 x 2 + a !3 x 3 + . . . + aij xj +.„+ ai„ x n = bi 
a 2 i xi + a 22 x 2 + a 23 X3 + . . . + a 2 j xj +...+ a 2n x„ = b 2 



aiixi + a i2 x 2 + a i3 x 3 + ...+ ayxj +...+ ai n x n = bj 
10 

ami xi + a^ x 2 + x 3 + . . . + a^ xj +...+ a™ x n - b m 

In the equations, "ay" is the attenuation coefficient of element u f at energy band "1"; 
"Xj" is the unknown projected content of element "/*; and V is the measured total attenuation 
at energy band V Variable "ay" is related to the total cross section over the source neutron 

1 5 energy range, which can be determined in advance. There are m equations (representing 

respective images at m energy bands) with n unknowns (corresponding to n elements that are 
mapped). When m > n , there exists a least squares solution. In order to reduce the number of 
equations (images), all elements other than H, C, N, and O are represented with one 
attenuation component, which is flat over the energy range. In the solution, the content of 

20 most other elements is split between "H" and "other." Therefore, larger relative error is 

expected for hydrogen, especially when a substantial quantity of aluminum is present. The flat 
component also compensates for internal scattering effects within the object. The linear 
attenuation model assumes the following two facts: (1) the attenuation by one element is 
linearly dependent on its projected content and (2) the existence of one element does not affect 

25 the attenuation of other elements. 

It should be noted here that although a linear model is used for fitting the data in the 
above description, it is not the only method that can be used. For example, non-linear fitting 
methods may be used or other fitting techniques, such as simulated annealing methods, may be 
used. A common aspect of all of these methods is that the entire spectrum of neutrons can be 

30 used for fitting and thus the utilization of neutrons is very efficient, thereby reducing both 
accelerator power and the amount of shielding used in a deployed system. 
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The equations, above, can be written in the compact matrix form, Ax = k The 
attenuation coefficients can be determined by experiment or by simulation. As long as m> n, 
the equation set is solvable. 

Consider the problem of finding the elemental composition vector xe R n , where R is 
5 the vector space, such that Ax- b when the attenuation coefficient matrix A e R"™ and the 
measured attenuation vector be R m are given. Usually, there are more equations than 
unknowns; that is, m>n 9 and the set is over-determined. Unless b belongs to R(A) t the vector 
space defined by A, the over-determined set has no exact solution. There are many possible 
ways of defining the "best" solution; one example is a least squares approach. 
10 The least squares solution of over-determined equations Ax - b is expressed as: 

A T A x = A T b xls - (A T A/ J A T b = a\ 

This does not mean that jc^ should be calculated in the above manner. In fact, this 
calculation is not recommended for two reasons: (1) calculating the inverse matrix is 
computationally intensive, and (2) it introduces large errors. Direct Householder transform 
15 can be used to solve the least squares problem without increasing the condition number. 

Householder transform reflects one vector on another vector, as shown in FIG. 4. The 
image vector, y, is determined by the object vector, x, and the rnirror vector, u. On the other 
hand, if we have the object vector, x, and want to reflect it into a certain direction (for 
example, y « || x || ei, where ei = (1, 0, 0, 0) T ), the mirror vector, u, can be derived from 
20 the reflection equation, y - x -2 (x«u) u - 1| x || ei . 

= x- H x H e t 
II He,) || 

That is, u can be found by subtracting from x the desired image, ||x|| ei. Normalization is then 
performed. In matrix form, 

y = (I-2uu T )x=||x|| e t 

25 where u is found in the above expression. The matrix, Q = I- 2 u u T , is called a Householder 
Matrix. It is orthogonal and normalized, and it has the following properties: 

Now it is straightforward to up-triangularize matrix, A, column by column (treat each 
column as an object vector, one at a time), until it consists of n up-triangularized equations and 
30 ro-n residual equations. The n equations can then be solved by backward elimination. 

9 
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The process is expressed as: Q„ ...Q 2 Q,Ax = Q n Q n .,...Q 2 Q,b, or 
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Q = Qn Qn-i 2; is also orthogonal and normalized. The resulting matrix, QA, is 
no more ill conditioned than matrix, A, k(QA)= k(A). In other words, the Householder 
5 transform does not increase the upper limit of error. Therefore, this method is, in general, 
expected to be more stable than other methods. 

In the detailed matrix expression above, * marks represent modified numbers and A 
represents small residual errors. The Householder method leads to the same solution as the 
normal equation does, x H = x T ls- When programming, there is no need to calculate the 
1 0 matrices, Ql , Q2, etc, explicitly. 

4. Setup and Neutron Source 

A typical fast neutron radiography setup is illustrated in FIG. 5. As noted, above, the 
neutron source is usually accelerator based. A proton beam or deuteron beam 32 is 
accelerated to the desired energy and bombards the target material 30 to produce fast neutrons 

15 by a nuclear reaction. Scintillation screens 68 are typically used for fast neutron radiography, 
and a mirror 70, a lens 72 and a CCD (chip) 74 camera are employed to record the image of 
objects 24. The mirror 70 is used to avoid direct neutron irradiation of the CCD camera 74, 
and the lens focuses light onto the CCD camera 74 and helps to collect more light. The light 
from neutron detection is extremely weak, so the scintillator 68, mirror, lens 72 and CCD 

20 camera 74 are light-sealed in a dark chamber 76. Additional shielding 78 is provided to 
protect the CCD camera 74 against direct irradiation from the neutron source. 

The most common way of producing the fast neutrons is through nuclear reactions, 
which can be triggered, as in the above-described apparatus, with a particle accelerator and a 
target. In particular examples, a p- 7 Li source or a deuteron-deuterium (DD) source can be 

25 used to produce the fast neutrons. The former source involves two nuclear reactions, while the 
latter source has only one reaction channel. 

p- 7 Li source: 7 Li + p^Be + n -1 .646 MeV 
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7 Zz + p^Be* + * - 2.076 M?K 

DD source: 2 /f + </= 3 /fe + n + 3.27 MeV 

The p- 7 Li reaction is endothermic. The resulting 7 Be nucleus can be in a ground state 
or in an excited state, emitting a 428 keV gamma ray when de-exciting. The threshold proton 
5 energy is 1 .88 MeV for the ground state reaction and 2.37 MeV for the excited state reaction. 
The cross section of the second reaction is much lower than for the first reaction when the 
proton energy is below 5 MeV, so the first group of neutrons can be thought of as mono- 
energetic contaminated by lower energy neutrons. 

The DD neutron source is the source most commonly used for neutron production in 

10 the 3 MeV to 12 MeV energy range. The reaction in a DD source is characterized by a single 
neutron group. The reaction is exothermic, and it has a Q-value of +3.27 MeV. The minimum 
neutron energy in the forward direction is 2.6 MeV. 

Neutrons of different energies are obtained by rotating the object-detector assembly 
around the DD neutron source, as shown in FIG. 6 and described further, below. It is 

15 estimated that a practical neutron flux can be achieved with IOOjiA average deuteron current. 

Using the angular dependence of the DD neutron energy spectrum allows the use of a 
fixed-energy accelerator, such as a radio frequency quadmpole (RFQ), and a thick target to 
generate fast neutrons. RFQ ion accelerators are compact, high-beam current devices usually 
designed to accelerate a specific ion to a fixed energy. Deuterons are the accelerated ions in 

20 one embodiment, and the target 30 can be high-pressure deuterium gas or a solid (deuterated) 
target. Solid targets can have a thickness of about 10 microns. A gas target in the form of a 
gas cell at 1-3 atm pressure can have a thickness on the order of a centimeter, though the 
optimal thickness is a function of the mass density of the target composition. In any case, the 
thickness of the target can be selected to produce a suitably broad range of neutron energies 

25 across the neutron-emission spectrum since the energy of emitted neutrons is a function of the 
energy of the accelerated deuteron, which loses energy as it passes through the target 
Consequently, the use of thicker targets results in a broader range of emitted-neutron energies. 
Where testing is focused on single peaks in procedures measuring attenuation by subtracting a 
detector reading at an off-resonance energy with an on-resonance energy peak, thin targets that 

30 reduce the energy of an accelerated deuteron by about 50 keV as the deuteron passes through 
the target are used. Whereas, in the broad-spectrum methods disclosed herein, the target can 
have sufficient thickness to reduce the energy of the accelerated deuteron by about 0.4 to 
about 0.8 MeV as the deuteron passes through the target. 
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It is advantageous to employ solid-state targets because they can be directly mounted at 
the end of a vacuum beam path, with deuterated material bombarded by a deuteron beam. 
Isotopes of hydrogen are readily absorbed in titanium or zirconium to form metal hydrides. 
These metals are usually evaporated onto a metal backing. Titanium is the most commonly 
5 used occluding metal because it has the lowest stopping cross section of elements with a 
strong affinity for hydrogen and a high thermal conductivity. A typical ratio of occluded 
deuteron atoms to titanium atoms is 1 .5, and the evaporated titanium layer ranges from 0.2 to 4 
mg/cm 2 . 

After passing from the target 30 through the object 24, neutrons are measured by a 

10 detector, which can include a scintillator 68 and a charge coupled device 74, as shown in FIG. 
5. In the scintillator 68, fast neutrons are scattered by hydrogen nuclei and the energetic recoil 
proton ionizes and excites the molecule or crystal of the scintillator, causing the emission of 
light. Suitable scintillators include plastic scintillators such as NE102 scintillators from 
Nuclear Enterprises, Inc. and BC400 and BC430 scintillators from Bicron Corporation 

15 (Newbury, Ohio, USA). Other suitable plastic scintillators are made by Applied Scintillation 
Technologies, Elgen and Scionix. Alternatively, the scintillator can be an inorganic ZnS(Ag) 
(silver-activated ZnS dispersed in polypropylene) scintillator. 

Scintillation light from the neutron detector is generally reflected with a mirror 70 and 
then collected with a lens 72 onto a charge coupled device (CCD) chip 74, which converts 

20 individual photons into electrons. A CCD is an analog integrated circuit that converts an 
optical image into an electronic one. Electrons, induced by photons are accumulated in each 
pixel during the exposure time and transferred (usually sequentially) to readout electronics 
thereafter. After a desired exposure time, the number of electrons in each pixel is converted to 
a digital level that is communicated to a computer coupled therewith, the computer generating 

25 the final digital image. 

In fast neutron radiography, the detector works in an integrating mode, in which the 
CCD camera 74 records the total amount of light emitted from small areas of the scintillation 
screen 68. Light output from individual neutron events is different The incident neutron may 
transfer its full energy to the recoil proton in a head-on collision or part of its energy when 

30 scattering at an angle. It could also scatter on a carbon nucleus, losing some energy, or it 

could experience multiple scattering within the scintillator. Furthermore, different amounts of 
light are expected for the same amount of total energy if the scattering histories (Le. 9 the 
number of prior collisions) are different Nevertheless, the sum total of emitted light at each 
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pixel relative to that at other pixels generally reflects the relative number of neutrons passing 
through the object at that pixel at a given neutron energy. 

The equation governing fast neutron radiography is km the exponential attenuation 

law. 

N = N 0 

In the above equation, "N" is the attenuated neutron number in a small area of the 
detector, "N 0 " is the corresponding source neutron number without attenuation, "n" is the 
effective attenuation coefficient (unit: cm 2 ) and V* is the effective total projected atomic 
content (unit: #/cm 2 ). 

Further description of neutron imaging conducted with a scintillator and a cooled CCD 
camera is provided in R.C. Lanza, S. Shi, and E.W. McFarland, "A Cooled CCD Based 
Neutron Imaging System for Low Fiuence Neutron Sources," IEEE Trans. Nuc. ScL, Vol. 43, 
1347 -1351 (1996), which is incorporated herein by reference in its entirety. 

The imaged object can be exposed to neutrons having different energies by rotating the 
aligned object-detector assembly around the neutron source, as illustrated in the embodiment 
of FIG. 6. A mount 22 supports the object 24 and positions the object 24 in the path of 
emitted neutrons. The mount 22 is coupled with a drive mechanism 28 (typically, a motor 
designed to generate physical displacement via circumferential rotation or other means) for 
repositioning the mount in various radial sectors of the total neutron-emission spectrum about 
the neutron-emitting target 30. Alternatively, the target 30 can be repositioned by the drive 
mechanism so as to put the object and detector at a different position within the neutron- 
emission spectrum (in such embodiments, the neutron-emission spectrum is rotated rather than 
the object and detector). The neutron-emission spectrum shown across the range of angles 
illustrated in FIG. 6 is a product of the use of a thick target 30, with the accelerated deuterons 
projected along the path illustrated by arrow 32. In one example, the object 24 and detector 26 
(which can comprise the components within dark chamber 76 in FIG. 5) are rotated by 
positioning the object and detector on a rotatahle wheel that is centered at the source and 
coupled with a rotary motor at its axis of rotation. 

This practice is different from tomography in which the detector and source are aligned 
and the object rotates around a fixed axis. Suppose the energy of the accelerated deuteron (Ed) 
= 2.3 MeV and AEd = 0.8 MeV (thick target). The energy range of emitted neutrons changes 
from 4.70 - 5.55 MeV at 0° to 2.27 - 2.33 MeV at 1 15°. FIG. 3 illustrates the attenuation, or 
neutron cross sections, for hydrogen 1 1, carbon 13, nitrogen 15 and oxygen 17 over a range of 
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neutron energies from 2 to 6 MeV. The illustrated boxes above the charted attenuation curves 
represent "energy bins" showing the range of neutron energies in this example at various 
angles shown in FIG. 6. At each nominal angle, the neutron energy spectrum varies over the 
image, so different attenuation coefficients are used for different parts in the images. For a 
5 measurement of hydrogen, carbon, nitrogen, oxygen and "other," neutron energy should 

ordinarily vary from as high as 4 .75 MeV to 5 MeV to cover the broad valley of nitrogen 15 to 
as low as 2.25 MeV to 2.5 MeV to cover the valley of oxygen 1 7. 

In mapping the concentration and distribution of two elements in an object, the object 
24 and detector 26 are positioned at at least three different angular positions within the 

10 neutron-emission spectrum, as shown in FIG. 6; and neutrons at respective energies 

corresponding with each of the angular positions pass through the object 24 to the detector 26, 
which generates a reading of the spatial distribution of the detected neutrons. Where, for 
example, the spatial concentration and distribution of each of hydrogen, carbon, nitrogen and 
oxygen in an object are mapped, detector readings at at least five energies are generated by 

1 5 exposing the object and detector to neutrons in the neutron-emission spectrum at at least five 
angular positions. The linear-attenuation model describe in section (3), above, can then be 
used to solve for the content of each element at each pixel within a larger two-dimensional 
image. 

With higher attenuation, fewer neutrons reach the detector array. The detector 
20 generates a reading of the number and distribution of neutron impacts upon the detector. This 
reading can be in the form of pixel-by-pixel measurements from a CCD camera, and it is 
communicated to a computer processor that executes software code stored in computer- 
readable memory. The software code performs the calculations, described above, for mapping 
elements based on input data from the detector readings. The processor, operating the 
25 software code, then generates a map or a plurality of maps of elemental compositions within 
the object and communicates that map to a display or printer so that the map can be viewed by 
the system operator. Alternatively, or in addition, the compositions and/or shapes of a 
substance or object to be detected can be stored in the computer-readable memory. For 
example, a library can be established for the compositions of various drugs or explosives. The 
30 processor can compare these compositions and/or shapes with the compositions and/or shapes 
derived from the detector readings to determine (if there is a substantial match) whether the 
substance or object to be detected is suspected to be present, whereupon additional attention 
can be devoted to the investigation of the suspect substance/object. 
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CCD defects lead to dark or white points, lines or clusters in the digital image. When 
radiation (neutrons, gamma rays or cosmic rays) hits the CCD chip, there might be an 
interaction with silicon, creating an energetic charged particle. The resulting charged particle 
will generate a lot of electron-hole pairs and leads to a white point or points, often referred to 
5 as "sparkles." The computer removes these sparkles with a nonlinear filter, which changes the 
white pixel to the median value of the pixels within its neighborhood the 3 x 3 or 5 x 5 
block surrounding it). Sparkles are detected with a deviation check, wherein a pixel value is 
replaced with the median value only when it deviates from the median within the 
neighborhood to a certain extent (e.g. t one standard deviation of the neighborhood values). 
1 0 This modified median filter removes sparkles but preserves fine structures in the digital image. 

SIMULATED EXEMPLIFICATION 

The following example is a simulation of elemental mapping using fast neutron 
resonance radiography in accordance with the methods described, above. 

In this simulation, the attenuation coefficients of the simulated object are generally 

1 5 smaller than the average cross section of the simulated object over the corresponding neutron 
energy range. Rather than calculating these coefficients directly from cross-section data, 
MONTE CARLO N-PARTICLE (MCNP) transport code was executed by a computer to 
simulate the path of neutrons into and through the object and the consequent attenuation of 
neutrons in the object, as explained in Judith F. Briesmeister, "MCNP — A General Monte 

20 Carlo N-Particie Transport Code," Los Alamos National Laboratory, March 1 997, which is 
incorporated herein by reference in its entirety. In the simulation for determining the 
attenuation coefficients, a slab of pseudo-single-element material was placed between a source 
and a detector, with the same geometry as shown in FIG. 6. The attenuation coefficients were 
calibrated for 0 - 1 .0 * 10 24 atoms/cm 2 for each element. 

25 Radiographic images were then simulated using the radiation transport code, "COG," 

developed at Lawrence Livermore National Laboratory, as explained in T. Wilcox and E. 
Lent, "COG - A Particle Transport Code Designed to Solve the Boltzman Equation for Deep- 
Penetration (Shielding) Problems," Vol. 1 , Users Manual, Lawrence Livermore National 
Laboratory Rept # M-221-1 (1989), and in J. Hall, "Monte Carlo Modeling of Neutron and 

30 Gamma-Ray Imaging Systems," SPIE, Vol. 2867, 465 (1996), both of which are incorporated 
herein by reference in their entirety. The source-to-object and object-to-image distances used 
in the simulations were 150 cm and 50 cm, respectively, giving an imaging angle of ± 10° and 
a magnification factor of 1,33. The energy spectrum of the neutron source was adjusted 
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downward from 5.55 - 4.70 MeV to 2.33 - 2.27 MeV in a series of 10 runs to simulate 
scanning the object through angles ranging from 0° to 1 1 5° relative to the axis of a real DD 
neutron source. Images of the direct (unscattered) neutron flux from the source and scattered 
neutrons in several bins at lower energies were generated using 200 x 200 pixel imaging 
5 detectors with appropriate energy and scattered-particle masks. These direct/scattered 

neutrons were summed together (weighted by energy) to form real images. The simulations 
were each run with 5 x 1 0 8 test particles on the ASCI Blue Pacific supercomputer at Lawrence 
Livennore National Laboratory using 256 processors operating in parallel. 

The object modeled in these simulations was a sort of 'terrorist overnight bag." The 

1 0 bag itself consisted of a thin aluminum shell (40 x 30 x 1 0 cm) with a wood handle, thick 
cloth covering and steel fittings. FIG. 7 shows a neutron radiographic simulated image of the 
bag at 0° (the 0° angle illustrated in FIG. 6) and marks the locations of hidden explosives and 
drugs. As shown in FIG. 7, the bag contained a newspaper, a bag of sugar (105g) 42, a stash 
of cocaine-HCl (105g) 44, a travel umbrella 46, a 10-cm switchblade knife 48, a paperback 

15 book 50, a block of plastic explosive (270g; 50/50 weight % mix of RDX and PETN), a pen 
and pencil set, a small camera 54, an automatic pistol 56 with extra ammo clip 58, a flat paper 
notebook and a selection of cotton, wool and nylon clothing items. The bag was heavily 
loaded and had an average density of around 0.5g/cm 3 . These items were not directly evident 
from the original images produced from simulated detector readings. 

20 The calculated elemental mappings of this simulation, derived using the attenuation 

model described above, are provided in FIGS. 8A-8F. In each of the simulated images, the 
elemental concentration was determined pixel by pixel as a function of the degree of neutron 
attenuation within the object at the region defined by each pixel. FIG. 8A illustrates the 
projected hydrogen content in the bag. FIG. 8B illustrates the projected carbon content in the 

25 bag. FIG. 8C illustrates the projected nitrogen content in the bag. FIG. 8D illustrates the 
projected oxygen content in the bag. FIG. 8E illustrates the summed projected content of all 
other elements in the bag. FIG. 8F illustrates the projected summed content of all elements in 
the bag. 

In the mappings, the simulated plastic explosive can be identified by its high 
30 nitrogen/oxygen content (see FIGS. 8C and 8D) and low hydrogen/caiton content (see FIGS. 
8A and 8B). The simulated stash of cocaine looks different from sugar in that the drug has 
about equal amounts of hydrogen (see FIG. 8A) and carbon (see FIG. 8B) but very little 
oxygen (see FIG. 8D). The bright bar overlapped on the drug in the carbon image (see FIG. 
8B) is the polystyrene handle of an umbrella. A pistol, a knife and a battery are visible in the 
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"other" picture (see FIG. 8E). A glass lens in the camera is also very clear in the oxygen 
picture (see FIG. 8D). Two aluminum buckles can be seen on the top in the hydrogen picture 
(see FIG. 8A). As we have mentioned, the calculation splits aluminum content into hydrogen 
and "other" results. 

5 

Although examples involving detection of explosives and drugs were emphasized, 
other materials can also be identified and mapped by the signatures in their elemental 
composition. Moreover, the NRR principles discussed herein are not restricted to the elements 
discussed in the examples (particularly, hydrogen, carbon, nitrogen and oxygen), but can be 

10 used with a variety of other elements, as well. Another area of interest toward which the 
methods and apparatus of this disclosure can be applied is in on-site mineral analysis. 

While this invention has been shown and described with references to particular 
embodiments thereof, those skilled in the art will understand that various changes in form and 
details may be made therein without departing from the scope of the invention, which is 

1 5 limited only by the following claims. 
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CLAIMS 

What is claimed is: 

1 . A method for mapping the elemental composition of an object, the method comprising 
5 the steps of: 

emitting a spectrum of neutrons from a source, wherein the energy of each 
emitted neutron is related to the angle at which the neutron is emitted from the source; 

positioning the object in the path of neutrons in a first radial sector of the 
neutron-emission spectrum; 
*0 detecting neutrons that pass through the object when the object is positioned in 

the first radial sector and generating a reading of the detected neutrons; 

positioning the object in the path of neutrons in a second radial sector of the 
neutron-emission spectrum; 

detecting neutrons that pass through the object when the object is positioned in 
1 5 the second radial sector and generating a reading of the detected neutrons; and 

measuring the presence of at least one element in the object as a function of the 
readings of neutrons detected in the first and second radial sectors. 

2. The method of claim 1, wherein the object is positioned in more than two radial sectors 
of the neutron-emission spectrum, and readings of detected neutrons are generated for 

20 each sector. 

3. The method of claim 2, wherein measurements for a plurality of elements in the object 
are generated from the readings of neutrons detected in the plurality of sectors. 

4. The method of claim 3, further comprising the step of evaluating the element 
measurements to determine whether an explosive material is present in the object. 

25 5 . The method of claim 3, further comprising the step of evaluating the element 
measurements to determine whether a drug is present in the object 

6. The method of claim 3, wherein the object is a luggage item. 
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7. The method of claim 1 , wherein the element is selected from hydrogen, carbon, 
nitrogen and oxygen. 

8. The method of claim 7, wherein the presence of a plurality of elements selected from 
hydrogen, carbon, nitrogen and oxygen is measured. 

9. The method of claim 1, wherein the object is positioned in different radial sectors by 
rotating the object about the source. 

10. The method of claim 1, wherein the source includes a target positioned relative to a 
particle accelerator such that neutrons are emitted from the target when accelerated 
particles interact with the target 

1 1 . The method of claim 1 0, wherein the target has a thickness sufficient to reduce the 
energy of an accelerated particle by at least about 0.4 MeV as the particle passes 
through the target. 

12. The method of claim 1, wherein the spatial distribution of the element(s) is mapped 
throughout the object. 

13. A method for mapping multiple elements in an object comprising the steps of: 

directing a first set of neutrons that have energies within a first range at an 
object, detecting the neutrons from the first set that pass through the object, and 
generating a reading of the detected neutrons from the first set; 

directing a second set of neutrons that have energies within a second range at 
the object, detecting the neutrons from the second set that pass through the object, and 
generating a reading of detected neutrons from the second set; 

directing a third set of neutrons having energies within a third range at the 
object, detecting neutrons from the third set that pass through the object, and 
generating a reading of detected neutrons from the third set; and 

generating a map of at least one element in the object from the readings of 
detected neutrons from each of the three sets. 

14. The method of claim 13, wherein maps of a plurality of elements in the object are 
generated from the readings of detected neutrons. 
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1 5 . The method of claim 1 3, wherein the steps for each set of neutrons are performed for 
more than three sets of neutrons, each set of neutrons having energies within a 
respective range. 

16. An apparatus for mapping elements in an object comprising: 

5 a particle accelerator for delivering accelerated particles; 

a neutron-emitting target positioned to interact with accelerated particles 
delivered from the particle accelerator to produce a neutron-emission spectrum; 
a mount for supporting the object to be imaged; 
a drive mechanism coupled with at least one of the mount, the particle 
10 accelerator and the neutron-emitting target to allow the mount to be at different 

positions within the neutron-emission spectrum; and 

a detector configured to measure neutrons that pass from the target through an 
object supported by the displaceable mount. 

17. The apparatus of claim 1 6, further comprising a computer processor in communication 
1 5 with the detector, the computer processor also being in communication with computer- 
readable storage medium storing software code for mapping a plurality of elements in 
the object based on readings from the detector. 

18. The apparatus of claim 16, wherein the target has a thickness sufficient to reduce the 
energy of an accelerated particle by at least about 0.4 MeV as the particle passes 

20 through the target. 
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